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Aba tract . During cooling, a hail of oxygen can fall in the center of a Carbon-Oxygen white dwarf. Conaequencea 
for white dwarfs evolution are considered. Formation of pulsars or type I Supernovae production can be the 
result of a difference in chemical composition . 

The following report is divided into four sections : (I) elementary information on the internal structure of 
white dwarfs, (2) basic information on highly correlated plasmas, (3) implications for phase separation in the 
core of cooling white dwarfs, and (4) consequences for evolution of white dwarfs. 

The present report can be considered as a continuation of a paper of Canal, Isern and Labay (1981). 

(0 Internal structure of white dwarfs . We just recall a few properties of white dwarfs which are relevant 
to the present problem. 

We describe a white dwarf as a non-rotating, spherical star, in hydrostatic equilibrium. The main problem 
concerns the equation of state. 


P- P( p, T ) 

which, in the major part of the star cannot be distinguished from the equation of state of a Fermi gas at 
zero temperature. 


P - P( P, 0) 

The region where the departure from T *^0 ^nfluen^gs the structure is extremely small, and represent a mass 

AM / M • 10*^’ T ' 
eri 

If we ignore the properties of the nuclei and introduce in the equations of hydrostatic equilibrium the 
general relativity corrections to the Newtonian theory of gravitation, we write 


dP C, \ I ( I ♦ ( P/ P c^)) (l-^CAffPr^/M c^)) 

P 

dr r^ (1-(2GM/ c^)) 

2 

where the term ( (1 (P/ pc)) represents the increase of inertial mass , the next term the increase in 

gravitationally attractive mass and the denominator term represents the change of the metric for Euclidian 
to non Euclidian . 

The knowledge of the equation of st^te is necessary in order to solve the equation of hydrostatic equilibrium. 
At high densities we can write the relativistic equation of state , 

P - KjC 

with \i^ ■ (A/Z) and 

- 1.23 lo'^ 

At lower densities, the non relativistic equation of state (neglecting electrostatic interactions) is 

P - K, ( p / 

1 2 

with K. j 9.91 . The change from one equation of state to the other takes place for a density of 

about to g cm 

For the kind of white dwarfs we are considering , the transition from the non relativistic to the relativistic 
equation of state takes place at a mass 

“rel -(4 n / G M^) ( p/ y K, « 8 . lo"^ 

M * 
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The mass radius relation for white dwarfs shows the existence of a maxinuD mass. For ■ 2 , we have the 
following properties of a white dwarf with the oaxinum loass : 

M - 1.36557 M (Sun) 

R ■ 1066 kn 

1.92 . 10*^ g cm ^ 

A white dwarf of high density , close to the limiting mass ( we should call it the Chandrasckhar-Kaplan 
limiting mass), has no other energy source than its internal thermal energy . As long as the inside of the 
white dwarf is fluid, the internal thermal energy is the sum of the thermal energy of the ions and of the 
thermal energy of the electrons , 

U - ( p/ A ) ( (3/2) k T ♦ ( 3/4) ( Z kV/ x) ) 

where x is related to the density by 

p - ( 8 tt/ 3 ) ( m^^ / h^ ) ( A m^ / Z ) x^ 

The corresponding cooling time scale is given by the heat balance condition . 

Again , for the kind of stars which we are considering, close to the Chandrasekhar - Kaplan limit, we have 
for the characteristic 'ooling time 

^ . ,o 6-056 ( T / 10®) yc»r«. 

This is valid however only as long as the white dwarf has not reached the soHd state. If we consider the ratio 
of the electrostatic energy to the thermal energy, 

r • ( / a k T ) 

where 

a-(3Z/4it N ) 
e 

is the ionic radius, the solidification taKes place when 

r - 160 . 

During cooling, the first point where solidification takes place is in the center of the star. Two quantities 
are imporcant to consider : (1) the melting temperature inside the star ; (2) the Debye temperature inside 
the star. It turns out that the heat content of the electrons is much larger tnan the heat content of the 
ions . 

When solidification takes place, the latent heat of solidification is liberated, and the consequence is a 
slowing down of the cooling . For the kind of stars which we are considering, this is completely negligible . 

2. Highly correlated plasmas . 

We have just seen that in a cold plasma , crystallization takes place when the ratio of the Coulomb energy 
to the thermal energy is about 160. We need not to consider in the following the quantum effects in the solid 

Let us consider now the case of a two component plasma. Carbon and Oxygen, which are the most abundant pro- 
ducts of the late stages of stellar evolution . 

There is a conjecture by Stevenson (1979) that an eutectic can exist , bet*>een Carbon an Oxygen, with the 
following characteristics : 

Xq (in number) - 0.332 

Tg - 0.628 T^ ( melting of Carbo ) 

■ 0.389 T^ ( melting of Oxygen ) 

with a density difference (Pollock and Hanae , 1974) : 

p(polid) - p(liquid) , 

• 3. lO’^ 

P 

and with A (Oxygen) « 15. 99468 and A(Carbon) ■ l2. 00000 , we obtain : 
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Pg(8olld)- p (liquid) 


3.10“* - 3.24 . lO"^ X 


1 - ( X / 4 ) 


;Xj, - X < 0.668 


p.(8olid) - p(liquid) 6.52 . 10~* - 4.07 . 10‘*X 
X > 0.668 

( 1 - (X/4)) 

As a consequence, we see that the density of the solid, eith^** of pure carbon or of pure oxy(.e.i , is 
always larger than the density of the liquid. 


During cooling, solid particles of Carbon , or solid particles of oxygen (according to the concentration) 
will appear and fall in the gravity field * When the chemical composition of the eutectic is reached, 
the solid particles of Carbon will be appreciably denser than the solid particles of Oxygen. In fact, 
we have : 


(Ap 

/p) 

( Carbon to Eutectic) 

- 3.91 - 

, 19-^ 

(Ap 

/p) 

(Oxygen to Eutectic) 

- 0.84 . 

, 19-^ 


The argument of Stevenson is that, if the cooling is fast enough, the eutectic solidifies; if it is very 
slow, it will experience separation of Oxygen and Carbon . However, due to the large difference in the 
buoyaicy forces. Carbon will settle faster. 


3. Plrise separation in white dwarfs . 

We shill first consider the temperature excess of the condensing solid. Following the argument of 
Jeffreys (1918) , and reproduced by Mason (1953) , we assume the microscopic diffusion is the way by 
which the atoms join the growing crystal. We then have: 

p(*olid) K - b (1 - X) N (A/Z) ni„. 

dt e H 

Similarly , we write that the latent heat of crystallization , released on the crystal is carried away by 
thermal conduc .ivity . We write, X being the heat conductivity, 

L R ^ - X (Tg - T ) 

dt 


from which we derive the temperature excess : 

T T . P (»-X) p L 
S A 

L is of the order cf 3 k T per ion. The microscopic diffusion coefficient is of the order of 

“ ■ ^ ''th 

The thermal conduc ivitv is large, heat being carried away by the relativistic electrons ( Scha: ’.man. 

White Dwarfs) . 

10 —1 ft 

For p • 10 g cm , T ■ 10 *K, this gives 

Tg- T * 1000 ''K ’« T . 

Let us now consider the speed at whi^h the crystals can fall inside the white dwarf. We consider that the 
force on the crystal is given by the Stokes formula. We then have: 


-(C M / r^) (4w / 3 ) ( Pjiquij) - 6 8 Pg V a < dr/ dt ) 

where v is the kinematic viscosity coefficient . 

If we consider that the radius of the crystal , a , is ..ime dependent, 

( a^/ 2) - D (l-X) t 


we obtain, 


^solid ^liquid 

( dr/ dt) - - (lb TT / 27) G ( D/v ) (A-X) t 

P 
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we can write 


r « Tq exr ( - ( t/ t) ) 

and derive the time scale of the phase separation , 


x' - (27 / 8 o ( P/ ( . / G ( V/ D) (1-X)-' 


The viscosity is due to the transfer of momentum to the ions « Therefore, v is of the order of D 
With p - 2 • 10^ g cm , Ap / p - 3.9 • 10 , we obtain , 


With Ap/p == 0.84 . 10 * 


T - 6 sec 


T - 13 sec. 


We can consider that the cooling time is much longer than the crystallization time and the infalling time . 
Two cases have to *e considered: 

(1) < 0.668 : Oxygen is predominant • 

Oxygen crystallizes immediately, as soon as the eutectic curve is reached. Carbon being soluble in oxygen, 
the temperature and the oxygen concentration decrease until the chemical composition of the eutectic is 
reached. The liquid, with the eutectic composition, which is left behind rsoves up and mixes with the liquid 
above, making it more rich in carbon. Further separation takes place until the whole solid oxygen corer is 
surrounded by the eutectic. 

Due to the slow cooling, crystals of C and 0 will form, but we can expect Carbon to fall faster than Oxygen. 
The melting temperature will have a tendency to rise above the melting temperature of the eutectic, with 
a tendency of the mixture to solidify quickly. The most reasonable assumption is that the oxygen corer is 
surrounded by a solid eutectic. 

(2) > 0.668 : Carbon is predominant . 

We have a similar picture, but, instead of an oxygen core, we have a carbon core, surrounded by the solid 
eutectic . 

4. White dwarf evolutiv 

(a) thermal effects. 

The infall of solid Oxygen (or solid Carbon) makes available a gravitational '^nergy which is of Che same 
order of magnitude as the internal heat content of the star, and larger chan the latent heat %diich is made 
available during crystallization . 

(b) role of accretion . 

For a white dwarf, close to the Chandrasekhar - Kaplan ILsit, the capture of a small amount of matter can 
make the star unstable. Capture at an appreciable rate can Cake place when the white dwarf belongs to a 
binary system. The effect is a pressure increase in the core of the star. 

Two cases have to be considered. 

(i) the white dwarf has an oxygen core. In that case, electron capture on oxygen can take place. 

Yhe decrease in Che number of free electrons per unit of mass generates instability . The, the star contracts 
about a a speed which is sort of average between the free fall time and the characteristic time of electron 
capture, until it reaches dynamical instability . 

What is expected is the formation of a neutron star. Hass ejection taking place at a later ^tage of 
collapse can lead to the escape of the star from the binary system. This uld explain the large value of the 
space velocities of pulsars in the galary . 

(ii) The white dwarf has a carbon core . Xrt t case, dynamical instability would take place before 
electron capture . In fact the explosive reaction , 


%iould take place even before the star has begun to collapse. The present conjecture is that a thermal 
runaway would take place and the star would explode as a type 1 supernova. 


225 



Bibl iography 


The theory oi the internal structure of White dwarfs can be found in 

S* Chandrasekhar : Stellar Structure, Chicago University Press, 1939 

There are several luonographies and review articles : 

E.Schatzman , White Dwarfs, North Holland Publishing Comp. 1958 
E.Schatzman , in Handbuch d.Phys. , Vol LI , 1958 , p.723 

L.Mestel , in Stars and Stellar Systems, Stellar Structure, L.H.Aller and D.B.Mc Laughlin Ed., Chicago 
University Press, 1965, p.297, 

A i-cent review paper is : 

H.Van Horn , Physics Today 19/9. 

Recent develoi;7»ents on the theory of highly correlated plasmas can be found in ; 

Strongly coupled plasmas, G. and Cavini P. Ed. Plenum Press, N.Y. and London 1977 . See especially the 

papers by 
Ichimaru S. p.189 
Schatzman E. p.407 

A colloquium , devoted to the problems of dense matter, includes many aspects of the properties of white 
dwarfs matter : 

Physique de la matiere dense , Journal de Physique , Tome 41 , C2 , 1980 . See especially the papers by 
D.J. Stevenson , p.61 
H.Van Horn , p.97. 

The conductivity of degenerate matter has been studied by several authors . See especially : 

A.Kovetz and G.Shaviv, Astronomy and Astrophys. 28, 315 , 1973. 

Flowers E. and Itoh N. , Astrophys. J. 206 , 218, 19“6. 

The viscosity due to ions has been estimated by 

H.Van Horn , in Low Luminosity stars, Kumar S, Ed, Gordon and Breach , N.Y. , London and Paris 1968, p.297. 

Concerning the idea of non-explosive collapse, see : 

R. Canal and E.Schatzman , Astronomy and Astrophys. 46, 229, 1976. 

The colloquium : 

White dwarfs and degenerate stars, l.A.U. Colloquium N'" 53 , H.M.Van Horo and V.Weidemann Ed. See especially 
R. Canal and J.Isern , p 52. 

R.Mochkovitch , These de troisieme Cycle, 1980 (tc be published ) 

The question of the effect of accretirn on white dwarfs has been considered again during the l.A.U. 

Symposium N® 93 : 

Fundamental problems in the theory *jf stellar evolution , D.Sugiraoto, D.Q.Lamb and D.N. Schramm Ed. Reidel 1981 


226 



